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Abstract: Tracy and Widom have evaluated the cumulative distribution of the largest eigenvalue for the finite 
and scaled infinite GUE in terms of a PIV and PII transcendent respectively. We generalise these results to 
the evaluation of EN{>^;a) (^U.iLixfloo,x]i^ - where X('2oo,a] = ^ ^ i-oo,X] and X('ioo,A] = ^ 

otherwise, and the average is with respect to the joint eigenvalue distribution of the GUE, as well as to the 
evaluation of ^Ar(A; a) :— YiiLii^ ~ Of particular interest are £'Ar(A; 2) and Fn^X; 2), and their scaled 

limits, which give the distribution of the largest eigenvalue and the density respectively. Our results are obtained 
by applying the Okamoto r-function theory of PIV and PII, for which we give a self contained presentation 
based on the recent work of Noumi and Yamada. We point out that the same approach can be used to study 
the quantities En{X] a) and Fn{X; a) for the other classical matrix ensembles. 
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1. Introduction and Summary 

Hermitian random matrices X with a unitary symmetry are defined so that the joint distribution of the indepen- 
dent elements P{X) is unchanged by the similarity transformation X i— > W XU for U unitary. For example, an 
ensemble of matrices with P{X) := exp^°^Q ajTT{X^) —: Hjli di^j) fo^' general g{x) > possesses a unitary 
symmetry. Such ensembles have the property that the corresponding eigenvalue probability density function 
p(Ai, . . . , Xn) is given by the explicit functional form 

1 ^ 

p{Xi,...,XN)^-^l[g{Xi) n l^fc-^jl'' (1-1) 

1=1 l<j<k<N 

C denoting the normalization. (Throughout the symbol C will be used to denote some constant, i.e. a quantity 
independent of the primary variable(s) of the equation.) The choice g{x) = , which is realized by choosing 
each diagonal element of X independently from the normal distribution N[0, l/\/2], and each off diagonal inde- 
pendently with distribution N[0, 1/2] + iN[0, 1/2], is referred to as the Gaussian Unitary Ensemble (CUE) and 
is the main focus of the present article. Specifically our interest is in the distribution of the largest eigenvalue, 
and the average values of powers (integer and fractional) of the characteristic polynomial Jli^i('^ ~ ^i) ^'^^ such 
matrices. 

Let i?Ar(0; (s, oo)) denote the probability that there are no eigenvalues in the interval (s, oo) for N x N GUE 
matrices. The distribution of the largest eigenvalue Pmax(5) is given in terms of En{0; (s,oo)) by 

Pmax(s) = ^^Jv(0; (s, oo)). (1.2) 
With i?Ar(s) specified by the solution of the nonlinear equation 

(i?^)2 -I- iiR'j^fiR'j^ + 2N) - 4{sR'j^ - Rn)^ = 0, (1.3) 



(an example of the Jimbo-Miwa-Okamoto a-form of the Painleve IV differential equation; see Eq. ( 2.18| ) below) 
subject to the boundary condition 



^^(^).:roo ^i/2(jv-i)! ' ^'-^^ 

it has been shown by Tracy and Widom |^ that 

£;Ar(0;(s,oo)) = exp(^- /" i?Ar(t)dt). (1.5) 
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The derivation in p5[ uses functional properties of Fredholm determinants (a subsequent derivation using the 
KP equations and Virasoro algebras has been given by Adler et al. In this work we will give a derivation of 



(1.5) based on the r- function theory of the Painleve IV equation due to Okamoto [g6|, and refined by Noumi and 
Yamada pl| . Our principal observation in employing this body of theory to problems in random matrix theory is 
that there is a deep and fundamental relationship between r-functions relating to the Hamiltonian formalism of 
the Painleve theory, and particular multiple integrals specifying averages with respect to the probability density 



function (1.1) in the case that g{x) takes the form of a classical weight function. From the random matrix 



perspective the classical weights are ||] 



Gaussian 



. x^-e ^ (a; > 0), Laguerre 
qix) = < , (1.6) 

^ {l-xY'{l + x)'' {-l<x<l),3acohi ^ 

(l + x2)-", Cauchy. 
To summarise the correspondence, which applies to any of the cases, we set out the following schematic table: 



Painleve Theory 


Random Matrix Ensembles 


r-fimction t[N]{s; a) 


Gap probability En{s) 
Averages En{s; a), Fjs[{s; a) 


Hamiltonian H[N]{s; a) 


Resolvent kernel i?7v(s;a) 
Logarithmic derivative of averages - 
UN{s;a),VN{s;a) 


Classical solutions - Weyl chamber walls 


Classical weights - Determinant structure 



Table 1. Correspondence between random matrix theory and Painleve theory. 



The quantity En{s] a) in Table 1 is specified by 

i?jv(A;a) := 



where X^'i^j^j = 1 for Aj G (-oo. A] and xf^ 



N 

n^('-co,A](A-A,r 

1=1 



(1.7) 



t(-oo A] ~ ^ ^"-"^ ^ "'^^^ A(_oo A] ~ ^ otherwise, and the average is with respect to the eigenvalue 

probability density function (|1 . Ij) . For general a we obtain the evaluation 



(1.8) 



En{s; a) — ENiso] a) exp i / UNit;a)dt 

J so 

(Eq. ( [4.14 ) with the substitution (4.10)), where UN(t;a) satisfies the nonlinear equation 

(U'j^f - 4(t[/^ - UNf + W^iJJ'^ ~ '2a)iU'j^ + 2N) = 0, 
(Eq. ( [4.15D ) subject to the boundary condition 

UN{t;a) 



t^-oo t 



1 



(1.9) 



(1.10) 
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(Eq. (4.18)). For (N+l) x (A^+1) dimensional GUE matrices pmax(s) is proportional to e En{s; 2). We therefore 
have 



Pmax(^) 



Pmax(so) 



exp( f [-2t + UN{t;2)] dtY 



(1.11) 



(Eq. O). 

The quantity Fpf^s; a) in Table 1 is specified by 

Fjv(A;a) := 



N 



1=1 



(1.12) 



(Eq. ( |l.7| )). For general positive integers a (1.12) has been computed by Brezin and Hikami in terms of the 
determinant of an a x a matrix involving Hermite polynomials. Note that for a not equal to a positive integer, 



(1.12) is well defined provided A has a non-zero imaginary part. For general a we obtain the evaluation 

i^w(A;a) =Fjv(Ao;a)exp(^ VN{t;a) dt^ , (1.13) 
(Eq. ( [4.33D ) where Vjv(i;a) also satisfies the nonlinear equation (4.15), but now with the boundary conditions 

(1.14) 



VN{t;a) ^ x^(l + 0(lA)) as t 



OO 



(Eq. ([4.35 )) where x = 1 for t — > oo and |x| = 1 for t — > — oo. In the case a = 2 this average is proportional to 
the polynomial part of the eigenvalue density for (iV+l) x (A^+l) dimensional GUE matrices, which in terms 
of the Hermite polynomial i?jv(A) is proportional to each of the 2x2 determinants termed Turanians [0 



ifjv(A) i/jv+i(A) 
H^(A) 



Hn+iW K+iW 
H'n+iW -H^7v+i(-^) 



Hn{X) Hn+iW 

Hn+i{X) i?Ar+2(A) 



(1.15) 



(which are of course proportional to each other) . The result ( 4.33| ) with a = 2 implies 



P(A) 



N^N+l 



Ni-^N+1 



exp 



-2t + VNit;2)] dt 



(1.16) 



(Eq. (p6|)). 

In Sect. 2 we review the r- function theory of the Painleve IV equation, revising relevant aspects of the work 
Noumi and Yamada 121 



of Okamoto 1 24,22], Noumi and Yamada ||2lj,^ and Kajiwara et al. |16|. The culmination of this theory from 
our perspective is the derivation of determinant formula expressions for the r-function corresponding to special 
values of the parameters in the Painleve IV equation. On the other hand, it follows easily from the definitions 
that Ef^f and can be written as determinants. These are presented in Sect. 4. The determinant formulas in 
fact precisely coincide with those occurring in Sect. 2, so consequently we can characterise both En and Fn in 
terms of solutions of the nonlinear equation ( 4.15| ). The theory presented in Sect. 2 also allows En, Fn to be 



characterised as solutions of a certain fourth order difference equation (Eq. (2.82)), and Un, Vn as solutions of 
a particular third order difference equation ( 4.19| ). 

Also of interest is the scaling Hmit of ^jLTj ) and ( |l.l2D with A i-^ V2N + X/V2N^/^. This choice of coordinate 
corresponds to shifting the origin to the edge of the leading order support of the eigenvalue density, then scaling 



Application of the T-Function Theory of PIV, PII to the GUE 



5 



the coordinate so as to make the spacings of order unity as TV ^ oo. We find the scaled quantities can be expressed 
in terms of particular solutions of the general Jimbo-Miwa-Okamoto a form of the Painleve II equation 







(Eq. ( ^.lOD ). Specifically, as already known from [ p6[ , 

E^°'\s) Jim £;^(0; V2iV + ^=^) = exp ( - ^ rit) dt) , 



N^oo 



(1.17) 



(1.18) 



where r{s) satisfies (5.10) with a = 0. Also 



lim (Ce-^'^/^ENis-^a)^ = E''°^'{sQ]a)e^p(^ f u{t;a)dt^, (1.19) 

■ , _ / /7:^T^ /fi -^Sn 



/2JV + s / \/2Af 1 / 



where u{s]a) satisfies (5.1C) subject to the boundary condition 

Aa^-l (4a2-l)(4a2-9) 



u{s; a) 



s^ — co 



1/4^^ + 



8s 64s4 
(Eq. ( ^.ll|) ). In the case a = 2 (Eq. (x7)) gives the formula 



+ 



Pmax(5) = Pmax(50) exp ( / ^(i; 2) dt 

J So 

(Eq. ( |S.20| )) for the scaled distribution of the largest eigenvalue in the GUE 
Analogous to the formula (5.7), for the scaled limit of Fn{^; a) we have 



(1.20) 



(1.21) 



F"°"(A;a):= lim 



/A 
v{t;a) dt^ , (1.22) 



(Eq. (5.31)) where f(s;a), like u{s;a), satisfies (5.10). The difference between u and v is in the boundary 
condition; for the latter we require 



v{t; a) 



-at 



1/2 _ ^ a(4a^ + l) 



(1.23) 



(Eq. (|3|)). The case a — 2 corresponds to the scaled eigenvalue density at the spectrum edge, which has the 
known evaluation ||l0| 



9^°f'(s) 



Ai(s) Ai'(s) 
Ai'(s) Ai"(s) 



where Ai(s) denotes the Airy function. In fact for all a G Z>o we have the determinantal form 



F^°f*(A;a) = (-l)"^''-!)/^ det 



d\3+^ 



Ai(A) 



j,/c=0,...,o-l 



(1.24) 



(1.25) 



(Eq. (p3|)). 

In Sect. 3 we present the r- function theory of the Painleve II equation in an analogous fashion to the theory 
presented in Sect. 2 for the Painleve IV equation. In particular we derive the second order second degree equation 
satisfied by the Hamiltonian (which is known from and j2^) as well as a fourth order difference equation 
satisfied by the t- functions. Also derived is the fact that the right-hand side of ( 5.33 ) corresponds to a r-function 



sequence in the PII theory, which is a result of Okamoto In Sect. 5 the results (q^), (^.7]), (5.31) and (5.33) 
are derived from a limiting process applied to the corresponding finite N results. 
A programme for further study is outlined in Sect. 6. 
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2. T-Function Theory for PIV 

2.1. Affine Weyl group symmetry. It has been demonstrated in the works of Okamoto p6t (in a series of papers 
treating all the Painleve equations), Noumi and Yamada (see also their works |2^,^,^) and the earlier 
work of Adler Q that the fourth Painleve equation 

y" = ^{yT + ly' + ^ty' + 2{e ~a)y + ^, (2.1) 
22/2 y 

can be recast in a way which reveals its symmetries in a particularly manifest and transparent form. 



Proposition 1 (|21,22|). The fourth Painleve equation is equivalent to the coupled set of autonomous difft 



er- 



ential equations (where ' — d/dt) 



fo = /o(/i - /2) + 2ao , 

/{ =/i(/2-/o) + 2ai , (2.2) 
/2=/2(/o-,/i) + 2a2 , 

with y — —fi and where the parameters aj G M with ao + ai + q;2 = 1 are related by 

a — ao — a2, P — —2a\ , (2.3) 
and the constraint taken conventionally as 

fo + fi+f2 = 2t. (2.4) 



Proof. Equation (2^) reduces the three first order equations of (2^) down to two. Eliminating a further variable 
by introducing a second derivative shows that y = —fi satisfies the PIV equation. The form of these equations 
implies 

(/o + /i + /2)' = 2ao + 2ai + 2a2 = fc, (2.5) 
fc 7^ constant, thus permitting the normalization given above. □ 

Note. Many differing conventions are in use for such a description of the PIV system and for example we have 
written 2aj (j = 0,1,2) in place of the aj used in j2^,^,^ in order to eliminate unnecessary factors of two 
appearing in the ensuing theory. 

The hyperplane ao + Q^i + 0^2 = 1 in parameter space {ao,ai,a2) G K'^ is associated with the simple roots 
ao,Q;i,a2 spanning the root system of type A'^K From this perspective the parameters ao,ai and a2 define a 
triangular lattice in the plane (see Fig. FI). 
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Let the fundamental reflections Si {i — 0,1,2) represent the automorphism of the lattice specified by a 
reflection with respect to the line Ofi = 0. Their action on the simple roots are given by 



(2.6) 



where the elements of the Cartan matrix 



A 



2 -1 -1 
-1 2 -1 
-1 -f 2 



(2.7) 



Let TT represent the lattice automorphism corresponding to a rotation by 120° degrees around the barycentre of 
the fundamental alcove C defined by a; > (i = 0, f , 2). Then 



(2.8) 



j e Z/3Z. The operators tt, Si obey the algebra 



Sj = 1, (SjSj + l)^ = 1, SjSj±lSj = Sj±lSjSj±l, TT^ = 1, TTSj = Sj+lTT 



(2.9) 



and generate W = (tt, sqj sij 52) defining an extension of the affine Weyl group associated with the ^2^-* root 
system. 



Proposition 2 (|20,23|). The Bdcklund transformations of the PIV system are given by the actions of the 



extended affine Weyl group W on the parameters as specified by (2.1) and \2.q) , and on the functions as specified 
by 



■^^{f,) = I, + 



2a, 



Hfj)=fj+i = 0,1,2) 



(2.10) 
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where the Uij are the elements of the orientation matrix 



U 



1 -1 
-1 1 

1 -1 



(2.11) 



associated with the boundary of the fundamental alcove \2l\ 



Proof. Let V denote one of n, sq, si, and let f3i := V{ai). Using and (2^) it's a simple exercise to 

explicitly verify that gi := V{fi) of the form ( p.lO[ ) satisfy the structurally identical equations 

g'o = 9o{gi - 92) + 2/3o , 

ff'i =ffi(ff2-5o) + 2/3i , (2.12) 
92 = 92(90 - 91) + 2^2 , 
thus giving rise to the stated Backlund transformation. □ 

Following p^jll^, in Tables 0, plthe actions (2.6), (2.8) and (2.10), are listed in tabular format. 





ao 


ai 


a2 


So 


-ao 


ao + ai 


ao + a2 


Sl 


ai + ao 


-ai 


ai + a2 


S2 


a2 + ao 


a2 + ai 


—a2 


TT 


ai 


a2 


ao 


Ti 


ao + 1 


ai — 1 


a2 


T2 


ao 


ai + 1 


a2 - 1 


T3 


ao - 1 


ai 


a2 + 1 



Table 2. Action of the generators of the extended affine Weyl group associated with the root system Aj^' on the simple roots. 





fo 


fi 


./2 


So 


fo 


2ao 

Jl + —r~ 

fo 


. 2ao 

J2 -T- 

Jo 


Sl 


, 2ai 

Jo 7~ 

Jl 


fi 


2ai 

J2 + —f- 

Jl 


S2 


2a2 
fo + —r- 

J2 


, 2a2 

Jl - 

J2 


/2 


TT 


/i 


/2 


fo 



Table 3. Backlund transformations for the PIV system. 
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From the earlier work of Okamoto it has been known that the PIV system, as for all the Painleve transcendents, 
admits a Hamiltonian formulation and that from this viewpoint the Backlund transformations are birational 
canonical transformations {q,p', H} i— > {q,p: H}. 



Proposition 3 (|26,16|). The PIV dynamical system is a Hamiltonian system {q,p;H} with the Hamiltonian 

H ^ {2p - q- 2t)pq - 2aip ~ a2q , 

(2.13) 

= 72/0/1/2 + "2/1 - ai/2 , 

and canonical variables q,p 



-fi = q, /2 = 2p . 
Proof. With H specified by ( 2.13| ), Hamilton's equations of motion read 



OH 

q' ^ ^ = q{^p~q- 2t) - 2ai 



dH 



p{2q -2p + 2t) + a2 



(2.14) 



(2.15) 



Substituting for p and q according to ( 2.14 ) shows that these equations are identical to the final two equations 



Note. Because — /i satisfies the PIV equation (2.1), it follows immediately from the first equation in (2.14) that 
q satisfies the PIV equation (p.l|). Furthermore, use of the first equation in ( 2.15| ) shows 



P^^[q' + q^ + 2tq + 2ai) , 



(2.16) 



so H is completely specified in terms of the Painleve IV transcendent ( |2.lD with parameters ( |2.3| ). 

There is a degree of ambiguity in constructing a Hamiltonian in that arbitrary functions of time can be added, 
and in fact there is a more symmetrical form 



Hs = 1/2/0/1/2 + y3(«i - «2)/o + 1/3(^1 + 2a2)/i - y3(2ai + a2)/2 



(2.17) 



which is central to the Okamoto theory (termed the auxiliary Hamiltonian). However for our purposes this 
complicates some later results so we prefer the unsymmetrical form. Furthermore, in the full theory of PIV | |2l| , 
^ the Hamiltonian Hq = H is associated with two additional Hamiltonians Hi ~ Tr{Ho),H2 = n'^{Hf)) but 
these are not required in the random matrix context. 

It is also true that H{t) can be specified as the solution of a certain second order second degree equation. 

Proposition 4 (|26,14|). The Hamiltonian { 2.1!\ ) satisfies the second order second degree differential equation 
of the Jimbo-Miwa- Okamoto a form for PIV, 



{H"Y - 4(iiJ' - Hf + m'(H' + 2ai){H' - 2^2) = 0. 



Proof. Making use of Hamilton's equations ( 2.15 ), we have for H{t) — Hit; q{t),p{t)), 

H' = /1/2 , 

H" = ./i/2(/2 - /i) + 2^2/1 + 2ai/2. 



(2.18) 



(2.19) 
(2.20) 
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Use of ( |2.13| ) and in ( ^.20D shows 



fi = 

f2 = 



-y^H" + tH' - H 
H' - 2oL2 ' 
i/aiJ" + tH' -H 



(2.21) 



H' + 2ai 

Substituting (2.21) in (2.13) gives the desired equation ( 2.18| ). □ 

For future reference we note that use of Tables ||, || shows that under the action of the generators of W ^ H 
transforms according to 



H 



2ao 



fo ' 

si{H) ^ H + 2ait , 
S2{H) = H -2a2t , 
7t{H) = H + f2- 2a2t. 



(2.22) 



2.2. Toda lattice equation. The r-function t = T{t) is defined in terms of the Hamiltonian H{t) by 

H{t)^j^\ogT{t). 



(2.23) 



It is possible to derive a Toda lattice equation for the sequences of r-functions {Tfe[n]}„=o,i,... {k = 1, 3) associated 
with the Hamiltonians 



H 



H 



(2.24) 



respectively (the reason for the subscripts 1 and 3 on r will become apparent subsequently). An essential point 
is that there exist shift operators from the algebra W which after n applications on H generate the shifts 
required by ( p. 24 ). There are in fact three fundamental shift operators [ pT| Ti := 7rs2Si, T2 := Si7rs2, Is := S2Si7r 
corresponding to translations on the root lattice by the fundamental weights Cjj, j = 1,2,3 of the root system 



^2^^ As can be checked from Tables 0, and (2.22) these operators have the property that 



TiH = H 



T^^H = H 



(2.25) 



Table ^ also shows that when acting on the parameters themselves, the same shifts occurring in the transformed 
Hamiltonian results, and thus 



Ti(ao, ai, Q!2) = (ao + 1, ai — 1, 02) 
T^^{ao,ai,a2) = (ao + 1,q;i,q;2 - 1) 



After a further n iterations the equations (2.25) can be written in the form 



-/(3)1 



(2.26) 



(2.27) 



where the subscripts (1) ((3)) refer to the system of Eqs. (2.2) with the parameters replaced as in the first 
(second) Hamiltonian (2.24) and use has been made of ( 2.13| ). 
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We remark that the two results of (2.25) are inter-related. Thus consider the mapping oj defined by multipli- 
cation by —1 together with the replacements 



(ao, ai,a2) ^ (-ao, -a2, -ai), (/o, /i, /2) ^ (-.fo, -.f2, -.fi)- 



(2.28) 



We see immediately that the system (2.2) is unchanged by as is the Hamiltonian (2.13), while we can check 
from Table that 



(2.29) 



Applying w to the first equation of (2.25) using (2.28) and (2.29) gives the second equation. With the r-functions 
Ti [n] and T3 [n] defined by 



application of (2.29) shows 



T^H=j^ log n [n] , T-"H = | log r, [n] , 



LUTi[n] = CT3[n]. 



(2.30) 



(2.31) 



In light of the relation ( 2.31 ), let us focus attention on the first equation of ( 2.25 ) only. 



Proposition 5 (|26,16|). The r-function sequence Ti[n\ corresponding to the parameter sequence {ao + n,ai 
n, a2) obeys the Toda lattice equation 



dt^ 



logai[n] 



ai[n+l]cri [n— 1] 



whe 



CTi 



[n] ■= Ce*'("i-"Vi[n] 



Proof. Following [£6,16 we make use of the first equation in ( 2.27 ) and consider the difference 

A crucial fact, which follows from Table |^ and (2.15), is that this difference is a total derivative 

d 



f{i)2[n]~T^ V(i)2N = ^log (/(i)i[?i]/(i)2[7i] -f 2(ai -71)). 



But it follows from (2.19) and (2.13) that 



/(i)iW/(i)2W + 2(ai - n) = l(rrff + 2t(ai - n)) = ^ log (e*'("-")ri [r 
and hence the right-hand side of ( 2.35| ) is equal to 

|,„,(il,og(e.-<".-...,[„,)). 

On the other hand ( ^.34 ) and (2.30) shows the left-hand side of (2.35) is equal to 



d_ / Ti[n+l]Ti[n-l] 
dt^^\ T^[n] 



(2.32) 



(2.33) 



(2.34) 



(2.35) 



(2.36) 



(2.37) 



(2.38) 
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Equating (2.37) and (2.38), and integrating shows that 



C 



ri[n+l]ri[ 



(2.39) 



and the stated result ( 2.32 ) follows. There is the ambiguity of a multiplicative constant C, possibly dependent 
on n but not on and this can be chosen freely, for example to render the Toda lattice equation in a simple 
form. □ 

The Toda lattice equation obeyed by the T3[n] with parameters (ao + ai, a2 — n) is obtained by applying 
the mapping uj to both sides of ( p. 39 ) and making use of ( p3l| ). This shows 



1 ( 



T-3[n+l]T3[71~l] 



which with 



gives the Toda lattice equation 



(73 [n+l]cr3[n-l] 



|jloga3[n] 



(T|[n] 



(2.40) 



(2.41) 



(2.42) 



Another way to deduce (2.40) is via (2.39) and the differential equation (2.18). Now, the second Hamiltonian 
in (2.24) is obtained from the first by simply interchanging ai and a2- On the other hand ai and a2 are 
interchanged in (2.18) if we replace t by it then replace H{it) by —iH(t). This tells us that in ( 2.39| ) we can 
make the replacements 



1 1-^ it, Ti{it) ^ 7-3 (i), ai ^ a2, 



(2.43) 



which indeed gives (2.4C). Furthermore, since (2.43) shows ti and T3 are simply related, it suffices to consider 
one sequence only, r^ln] say. 



2.3. Classical solutions. For a special initial choice of the parameters it is possible to choose T3[0] = 1, and then 
to determine T3[1] in terms of a classical function, that is to say the solution of a second order linear differential 
equation. What is essential here is the condition for the decoupling of the two independent first order differential 
equations so that what remains is a Riccati equation. 



Proposition 6 (|26|). For the special initial choice of the parameters such that = 0, i.e. for parameters 
(1 — ai, Q!i, 0) the first nontrivial member of the t -function sequence t^[V\ satisfies the Hermite- Weber equation, 



<[l]=.-2^[l]-2air3[l]. 
Proof. With n = 0, {1 ~ ai + n,ai, —n) implies 02 — 0. Now we see from ( ^.13 ) that 



H 



Q!2=0 



pq{2p — q — 2t) ~ 2aip 



(2.44) 



(2.45) 



which allows us to take 



H 



Ct2=0 



0, 



(2.46) 
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provided we set p = 0. Recalling (2.30) this implies 

T-3[0] = 1 . 

We read off from the n = 1 case of the second equation in ( ^.27 ), together with (2.14), that 



H 



and thus, after recalling the second equation in ( 2.30 ) and ( 2.46 ), 



^logr3[l]=g. 



(2.47) 



(2.48) 



(2.49) 



The first of Hamilton's equations ( 2.15| ) gives, with 0:2 = and p = (after the differentiation), the Riccati 
equation 



q' = -q^ - 2tq - 2ai. 



(2.50) 



Substituting (2.49) this reduces to the linear equation ( ^.44 ) first obtained in the present context by Okamoto 



Proposition 7. Two linearly independent solutions to the Toda lattice equation {2.4-L) for sequences of t- 
functions with parameters (ao + n,ai,—n), n > 0, starting from the Weyl chamber wall a2 ~ are given 
by the determinant forms 



T3[n](i;ai) = Cdet / {t - x)-"'+'+^ e-""' dx 



and 



f3[n](t;ai) =Cdet / {t - x)-°"+'+^ e-"" dx 

J —00 

Proof. In the special case ai — 0, we observe that a sohition of ( ^.44 ) is 



z J— 0,...,n — 1 



-i J— 0,...,n— 1 



r3[l] = C f 

J — O 



(2.51) 



(2.52) 



(2.53) 



In fact it is possible to solve ( 2.44 ) in a form analogous to ( 2.53| ) for general ai. Thus consider the integral 



Ia{t) := / {t-x^e-^ dx, 



(2.54) 



and suppose temporarily that Re(a) > — 1. Simple manipulation gives 

lait) = tla-l{t) + 1/2 /* it - xr-'^e-^" dx 



^tIa-l{t) + ^—^Ia-2{t). 



But 



(a - l)Ia-2{t) = -^Ia{t), la-lit) = - ^M^) 

a at'^ a dt 



(2.55) 



(2.56) 
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Thus we see that /_Q,j(i) satisfies (2.44) and this imphes 



T3[l] = C / {t-xy'e-'^ dx 



(2.57) 



where we require Re(Q;i) < 1. 

Starting with ( ^.47 ) and ( ^.57 ), up to a multiphcative constant the r-functions r^j^n] [n — 2,3,...) are 
uniquely specified by the Toda lattice equation (2.42). In fact it was known to Sylvester (|l^ pp. 115-117) that 
the solution of ( 2.42 ) with initial condition 173 [0] = 1 is the double Wronskian or Hankel determinant 



a^[n\ — det 



^3[1] 



(2.58) 



0,...,n — 1 



Recalling ( 2.41 ) and ( 2.57 ) we therefore have 
T3[n] T3,[n]{t; ai) Cdet 



; d^+3 
dP+3 



(2.59) 



i,j—0,...,n — l 



Makiiiff use of ( |I55| ) we can check that 

;(e*' f {t- a;)-"ie-^' dx^ = 2^6*' / {t - xy'+Pe-""' dx 



(P_ 

'dtp' 



(2.60) 



so (2.59) can also be written in the final form of (2.51). 

The second linearly dependent solution of (2.44) can also be written in an integral form similar to (2.54). 
Thus we see the integral 



{t-xfe-^'' dx, 



(2.61) 



satisfies the formulas (2.55) and (2.56), and thus satisfies (2.44) with a = —ai. Hence in addition to (2.57) we 
have the solution 



/oo 
{t-x)-"^e-''^ dx , 
-00 



(2.62) 



(note that for ai not equal to a non-positive integer, this is well defined only if i ^ R). Proceeding as in the 
derivation of ( ^.5l| ), we deduce from the Toda lattice equation (2.42), and the initial values (2.47), (2.62), the 
sequence of r-functions given by ( 2.52 ). □ 

Let us now consider the sequence of r-functions fi [n] . 



Proposition 8 (|26|). The sequence of t- function solutions to the Toda lattice equation (2.3^) fi[n], n>Q, 
corresponding to the parameter sequence (ag -I- — n, starting from the line ai = has the determinantal 
form 



fi[n]{t;-p) = Cdet 



p+i+j 



i.j—0.....7i — l 



(2.63) 



for -a2 = P e Z>o- 
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Proof. This sequence can be obtained from T3[n] by the (inverse of) the mappings (2.43). Replacing t by —it in 



(2.51) does not lead to an integral of interest in random matrix applications, but doing the same in (2.52) gives 



fi[n](i;a2) = Cdet / {t - ix)-°'^+'+^ e'^ dx 



i j'— 0,...,n — 1 



(2.64) 



which is of interest. We recall that for rj;"'' the parameters (ao,ai,a2) hi the corresponding Hamiltonian are 
given by 



(1 - q;2 + n, -n,a2). 



(2.65) 



For p G Z>o we know 



{t - ixfe-'' dx = y^2-PHp{t), 



(2.66) 



and thus setting a2 — ~p equation ( 2.64 ) yields ( 2.63| ). □ 

Note that with p — N , n — 2, this is precisely the final determinant in ( 1.151 ). 



2.4- Bdcklund transformations and discrete Painleve systems. It has been known that some of the Backlund 
transformations of the PIV transcendent can be identified with discrete Painleve equations |^,|llj, although no 
systematic study has been undertaken for this class. We will find that the difference equations for fj[n], H[n],T[n] 
which are generated by the Backlund transformations for the two shift operations Tj^^,Ti are in fact manifes- 
tations of discrete Painleve equations. 



Proposition 9 (|11|). The Bdcklund transformations of the PIV system corresponding to the shift operator 
T^^ generating the parameter sequence of (ao + n,ai,~n) with ri G Z, < ao,ai < 1 and ao + ai = 1 are 
second order difference equations of the first discrete Painleve equation dPI, namely 



Xk+i +Xk +Xk-i = 2t + 



(y2 + «i) + (-i)'(y2-«i) 

Xk 



k > 1 



(2.67) 



whe 



X2n+1 = /(3)2N, X2«+2 = /(3)oM ">0 



(2.68) 
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Proof. The action of the shift operators on the fj is expressible in a terminating continued ft'action, which for 
T3 and its inverse takes the form 

Tsih) = ,/i - ^ , (2.69) 

2(ai+a2) 

r3(/i) = /2 + - (2.70) 

h 

2a2 2(ai+a2) 

T,{h)^h + -, ^ (2.71) 



T^\h) = h-^+ — — ^ , (2.72) 
Jo 2ao 

fi + ^ 
Jo 

. 2(Q!o + ai) 

T^\.h) = fo- , (2.73) 

2ao 
h + —r~ 

JO 

2ao 

T^\h)^h + ^. (2.74) 
Jo 

as one can verify using the action of the afhne Weyl group reflections and diagram rotations as given in Tables 
H, ^. For simplicity of notation we suppress the subscript (3) labelling the sequence (ao + n,ai, —n) during the 
discussion of our proofs as there is no risk of confusion. Taking the first and last members of this set, now at 
the n*^ rung of the T3 ladder, and adding their unshifted /-variable we have 

2n 

fo[n] + h[n-l] ^2t- }2[n] + —- 

, (2.75) 

+ /2N = 2t- /oN + , 
Io[n\ 

so that one has a closed system. This can be recognised as the two components of a staggered system of difference 
equations and employing the definitions of Xk above we arrive at the discrete Painleve equation dPI. □ 

In terms of the coordinate and momenta of the Hamiltonian system this difference system was found by 
Okamoto [^|2^ and can be expressed as 

q[n + l] = {2t + q[n]~2p[n]) \ ^ i o n ' 2.76 

2(n+Q;o) - q[n\(2t+q[n\-2p[n\) 

q[n-l\ = -2p[n]^iUfLl , 2.77 

p[n + l] = -i/,q[n] + ^,^7 ^"9 [ 1 ' (^.78) 
' 2t + q[n\ — 2p[n\ 

^ ^ ^ q[n]p[n]-n ^^^^^ MM - a. (^^^g^ 
2p[n\ q[n\p[n\ — n 

Consequently a third order difference equation exists for the Hamiltonian through the relation 

H[n + 1] - H[n] = . (2.80) 



h 2a2 

h--r- 
h 

2ao 2(ao + ai) 
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Eliminating p[n\ between ( 2.7(^ ), (2.77) we find a second order difference equation for q[r 



nq[n]q[n-l](At + 2q[n] + q[n+l] + q[n-l] 



= 2 



{n + l)q[n + l] - nq[n-l] - (2t + q[n] + q[n+l])(ai + i/2q[n]{2t + q[n])) 
{n + l)q[n + l] - nq[n-l] - {2t + q[n] + q[n-l]){-ai + i-/2q[n]{2t + q[n] + q[n+l] + q[n-l])) 



(2.81) 



Use of — H[n] = q[n] leads to the third order equation in H[n]. 

In addition we have a higher order difference equation for the r-function. 

Proposition 10. The r-function sequence, appropriately normalised, associated with the shift operator T^^ with 
parameter values (ao + ai, —n), n G Z>o, < ao, ai < 1 and ao + ai = 1 satisfies the fourth order difference 
equation 



At^ {2nT'^[n] - T[n+l]r[n-l]^ x {2{n-ai)T'^[n] - r[n+l]r[n-l] 
X {T[n-2\T[n + l]T[n] + 2r^[n-l]r[n+l] + 4n(ai-n)T^[n]T[n-l] 
X {T[n+2]T[n-l]T[n] - 2r^[n+l]r[n-l] + 4n(ai-n)T^[n]T[n+l] 



= |r[n+2]r[n-2]r^[n] - 167i^(ai-n)V^[n] 

+ 16n(ai-n)(ai-2n)r^[n]r[n+l]r[n-l] - 4(2n2-2ain+l)r[n]T2[n+l]r^[n-l] 
+ r[n+2]r2[n-l] (^r[n+l]T[n-l] + 2{ai + l~2n)T'^[n]^ 

+ r [n - 2] + 1] (t [n + 1] T [?i - 1] + 2 (ai - 1 - 2n)r2 [n] 



(2.82) 



Proof. We first seek to express all the fundamental quantities in terms of the product /iM/2M- By multiplying 
the two transformations (2.69) and ( p. 70 ) we find a quadratic relation for fi[n] (and 1]), 



(2.83) 



(2.84) 



/i[n](2t - /i[n]) - /i[n+l]/2[n+l] + fMf^H + 2ai . 
Next we multiply ( ^.70 ) by fi[n] which yields a relation for the product 

/i [n] fi[n-l\^ fi [n] f2 [n] . 

/1N/2N + 2n 

One can verify then that a linear proportionality exists between fi[n\ and /i[n— 1] via the product /i [ri] /2 [jt^] , 

/iN I/1N/2N + + 2a^ - ^ [n] [n] 1 

I +271 J ^^^^^ 

= |/i[n + l]/2[n + l] + /1M/2W + 2a, - f,[n]f2[n]i^j^(^j^±^] , 

L /1N./2N + 2n J 

so that fi[n\ and /i[n— 1] may now be linearly related to fi[n\f2[n\. Multiplying these two later relations and 
using 



T3[n+l]r3[7i-l] , f r 1 f r 1 



(2.86) 



18 



P.J. Forrester, N.S. Witte 



with C = 1 to introduce the r-functions we arrive at ( |2.82| ). □ 

Note. The difference equations ( ^.81 ) and ( ^.82 ) have the advantage of being of the lowest order we have found 
possible, but the disadvantage of not being linear in the highest order terms {q[n + 1] and T[n + 2] respectively). 
In fact difference equations linear in the highest order terms can be given by increasing by one the order of the 



equations in each case |28 



Applying the operator lu (recall (2.2?)) we obtain analogous results for the sequence generated by Ti 



Proposition 11 (|11|). The Bdcklund transformations generated by the shift operator Ti corresponding to the 
parameter sequence {uq + n, — n, 02) with n G Z, < ao, 02 < 1, o-nd ao + 02 ~ 1, are second order difference 
equations of the first discrete Painleve equation dPI, that is 



Vk+i +rik + Vk-i = 2i 



k-[l + (-ma2-y2[l-(-m 



rik 



fc > 1 , 



whe 



V2n+1 = /(l)lN , V2n+2 /(l)oN 



n > 



(2.87) 



(2.88) 



Proof. This follows immediately upon applying w to both sides of (2.67). □ 

The analogue of (2.81) for the parameter sequence generated by the shift Ti can be found by applying the lo 
map to this relation, 

— 2np[n]p[n— 1] ^2t — 2p[n] — p[n + l] — p[n~l]j 

= {n + l)p[n + l]-np[n-l] + {t-p[n]-p[n + l]){a2 + 2p[n]{t-p[n])) (2.89) 

X {n+l)p[n+l] — np[n—l] + {t — p[n] — p[n—l]){—a2 + 2p[n] (t — p[n] — p[n + l] — p[n—l])) , 

and this implies a third order difference equation for the Hamiltonian via 

H[n+1] - H[n] = f2[n] = 2p[n] . (2.90) 

There is also a higher order difference equation for the r-function which can be derived using the relation 

Ti[n+l]Ti[n-l] 



fi[n]f2[n]-2n , 



(2.91) 



although we do not reproduce this here. 



3. T-Function Theory for PII 

3.1. Affine Weyl group symmetry. In the general Painleve theory the second Painleve equation naturally appears 
as a coalescence limit of PIV. From the work of [|6| it is known that in random matrix theory PII occurs in 
the edge scaling limit of the GUE. This suggests that before studying this limit we should develop a theory of 
PII analogous to that developed for PIV in the previous section. We take the PII equation to be defined in the 
standard manner 



y" ^2y^ + ty + a 



(3.1) 
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Proposition 12. The second Painleve equation with the transcendent y = q{t) and parameter a is equivalent to 
the system of first order differential equations 

f^ = ~2qfo + ao , (3.2) 
f[ = 2qfi + ai , (3.3) 

where fo + f\ = 2q^ + 1 and + ai = 1 with a = ai — 1/2 = 1/2 — ao- 

Proof. This is established by ehminating p through the substitutions fo — 2q^ — p + t and fi=p. □ 

Proposition 13. Let W = (sq, Si,7r) be the extended affine Weyl group of the root system of type A^^^ generated 
by the reflections So,Si and the diagram rotation tt, with action on the roots ao,ai as given in Table^, 





ao 


Oil 


So 




ai + 2ao 


Si 


ao + 2ai 


-ai 


TT 


ai 


ao 



Table 4. Action of the generators of the extended afiine Weyl group associated with the root system A^^' on the simple roots. 



The coupled system {3.i), (3.S) is symmetric under the Bdcklund transformations induced by the elements of 
the above affine Weyl group as specified in Table 





fo 


fi 


q 


So 


fo 


Aaoq 2al 
Jo Jo 


'^^ T 
Jo 


Si 


Aaiq 2a\ 
Jo+ , + ,2 
Ji fi 


fl 


fl 


TT 


fi 


fo 


-q 



Table 5. Backlund transformations for the PII system. 



Proof. This can be directly verified using the equations of motion (3.2), ( |3.3| ). □ 
Underlying the dynamics of the PII system is a Hamiltonian structure. 

Proposition 14. The PII dynamical system is equivalent to a Hamiltonian system {q,p]H} with Hamiltonian 

H = -1/2/0/1 - aiq , (3.4) 

and canonical coordinates and momenta q,p defined by 

p^fi, 2q^ = fo + fi-t. (3.5) 



20 



P.J. Forrester, N.S. Witte 



Proof. Using the symmetrised differential equations (B.2 3.3) the Hamilton equations of motion 



q' =p- - 1/2* , P' = "^qP + Oil , 



(3.6) 



can be verified. □ 

Remark. The fundamental domain or Weyl chamber for PII can be taken as the interval a S (— y2jO] or 
a G [0, 1/2), and so there exist identities relating the transcendents and related quantities at the endpoints of 
these intervals. In particular, denoting the transcendent q{t, a) and with = 1, t = we have |l| 



9(t,i/2e) - 62-1/3 1^5(s,0) . 

q(s, 0) as 

The action of the affine Weyl group on the Hamiltonian is given in Table ^, 



(3.7) 





Ho 


Hi=Ti{Ho)=Ho + q 


So 




Hi 


Si 


Ho 


Hi + ^ 
Jl 


TT 


Hi 


Ho 



Table 6. Backlund transformations of the Hamiltonian. 



We define two shift operators corresponding to translations by the fundamental weights of the afRne Weyl 
group A\^\ 

Ti = vrsi, T2 = siTT , (3.8) 
although T1T2 — 1 so only one is independent. Their action on the parameter space is given in Table ^ 





ao 


ai 




ao + 1 


ai — 1 


T2 


ao - 1 


ai + 1 



Table 7. Action of the shift operators on the simple roots of the root system A-[ 



(1) 
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Proposition 15. The Bdcklund transformations corresponding to the shifts are given by 

(3.9) 

(3.10) 
(3.11) 

(3.12) 

(3.13) 
(3.14) 

(3.15) 

(3.16) 

Proposition 16. The Hamiltonian H{t) satisfies the second order second degree differential equation of Jimbo- 
Miwa-Okamoto a form for PII, 

{H" f + 4 {H'f + 2H'[tH' -H]- y^aj = . (3.17) 
Proof. Using first two derivatives of H 

H' = -1/2/1 , 

(3.18) 

H" = -qfi - 1/20:1 , 



Tiifa) 


= fi~ 


f + 
Jo 


2a 

Tl 




— Jo , 






T2{fo) 


= fi , 






T2{fl) 


= /o + 


Aa^ q 
— — + 


2a 


h 


T? 


Ti{Ho) 


= Ho^ 


-q , 






= Hi - 


-q + Y' 

Ja 




T2{Ho) 




h 




T2{Hi) 


= Hi - 


- q ■ 





one can solve for q, /o, /i and tlien substitute these back into tlie expression for H. The result (3.17) then follows 
after simplification. □ 



With H given by (3^), and p and q specified by ( |3.q ), Hamilton's equation for q' implies 

p = q' + q^ + l/^t . (3.19) 
Thus H can be expressed in terms of the Painleve II transcendent q according to 

H - 1/2(9')' - + 'htf - (a + 1/2)9 . (3.20) 



3.2. Toda lattice equation. The r-functions are defined as before ( 2.23 ) and corresponding to each sequence 
generated by the shift operators is a Toda lattice equation. 

Proposition 17. The r-function sequence generated by the shift operator Ti with the parameter sequence (ao + 
n, ai — n) for n > 

T^{H) = H[7i] = j^\nT[n], (3.21) 



obeys the Toda lattice equation 



T[n+l]r[n-l] d% . , 
C — = 37^1n^N ■ 3.22 
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Proof. This parallels the argument employed for PIV case, by utilising the relations 



H[n + 1] - H[n] ^ q[n] , 



(3.23) 



and 



along with 



(3.24) 



(3.25) 



□ 



3.3. Classical solutions. When the parameter values are those on a chamber wall (a point) ai = n G Z then the 
T- functions are known to be expressible in terms of Airy functions |^ . 

Proposition 18. The solution for the first non-trivial member of the r-function sequence r[l](t) generated by 
the shift operator Ti with initial parameters (ao)Cki) = (I7O) that is bounded as t ^ —00 is 



(3.26) 



The n member of this sequence is 



T[n]{t) = Cdet 



(3.27) 



J i,i=0,... ,n-l 

Proof. Starting from ai = at n = one can take p[0] = so that H[Q] — and conventionally t[0] — 1. Using 



(3.13) we find that H[l] = q[0] and so the equation of motion (3.(:) gives the second order linear differential 
equation 



r"[l] + i/2ir[l] =0 



(3.28) 



and thus (3.26). The determinant formula (3.27) follows from (2.58). □ 



Another special parameter value of the Hamiltonian system (3.4) is ai = —1/2 |l37| when Hamilton's equation 



(3.6) permit the solution 



iq.p)^it-\i/2t) 



(3.29) 



However the corresponding value of H is not zero so in this case we do not have r[0] = 1 (rather H[0] and 
thus logr[0] is a rational function of t), and thus the sequence of T-functions generated by Ti is not given by a 
determinant Nonetheless the Backlund transformations of Prop. 15 show that H[n] and thus log T[n] remain 
rational functions of t for all n = 1,2,3,.... 
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3.4- Bdcklund transformations and discrete dPI. The discrete dynamical system generated by the Backlund 
transformations is also integrable and can be identified with a discrete Painleve system 

Proposition 19. The members of the sequence {q[nW, n > generated by the shift operator Ti with the param- 
eters (ttQ + n,ai — n) are related by a second order difference equation which is the alternate form of the first 
discrete Painleve equation, a-dPI, 



a + y. 



a — 1/ 



g[n]+g[n-l] + 
Proof We deduce from and Table | that 

a -1/2 



= -2q^[n] ~ t 



Ti{q) = -q- 
72 (g) = -q 



p-2q^-t 
P 



so eliminating p through the combination of these two we arrive at the stated result. □ 
The full set of forward and backward difference equations are [§61 



Q - y2 - " 

p[n\ — 2q[n\'^ — t 
a + 1/2 - n 



q[n+l] = -q[n] + 

q[n-l] = -q[n] - 

p[n\ 

p[n + l] = -p[n] + 2g[n]2 + t , 
p[n-l] = t-p[n] + 2 (q[n] + 



a + 1/2 - n 



p[n\ 

The discrete Painleve equation ( 3.30| ) implies a third order difference equation for the Hamiltonian 

+ 1/2 ~ a — 1/2 — n , , , , , , T 

„r ^ rr. + ' . = "2 if n+1 - H[n]f - t , 

H[n+1\ — H[n—1\ H[n+2\ — H[n\ 

because q[n] — H[n + 1] — H[n]. Equations ( |3.35D and ( 3.36D also imply 



1 a 

q[n] = p[n] {p[n-l\ - p[n+\]) - 

4a„ 2p[n\ 



a + 1/2 - n. 



Using this to eliminate q[n] (in say (3.35)) yields a second order difference equation for the p[n\ 

^.2 



-^pM {p[n+l]-p[7i-l]f + ^ 
Aaf^ p-^ [n\ 



2p[n]-p[n + l] -p[n-f] + 2t = 0. 



Furthermore, Eqs. ( |3.2l| ), ( |3.22D and (|3.25| ) give 



^r[n+^r[n-f]^^j^j 



So substituting in (3.39) (with say C = —2) implies a fourth order difference equation for T[n], 
^(r[n-2]r2[n+f] - T[n+2y[n-l]f + y^alT^[n] 

+ r [n - 2] r 3 [n] r 2 [?i + 1] + T [n + 2] [n] [71 - 1] + 2r 3 [n - f ] [n + f ] + t [n] r 2 [?7 - f ] r 2 [n + f ] = 



(3.30) 

(3.3f) 
(3.32) 



(3.33) 

(3.34) 
(3.35) 
(3.36) 

(3.37) 

(3.38) 

(3.39) 

(3.40) 



(3.4f) 
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While (3.39) and the corresponding equation for r[n] provide a polynomial relation between the smallest set 
of consecutive sequence members and {t[?t.]} we have found possible, they have the disadvantage of not 

being linear in the highest order term and T[n+2] respectively). This disadvantage can be remedied by 

increasing by one the order of the equation in each case. Thus by replacing n by n — 1 in ( 3.3j ), then adding 
the result to the original equation and using (3.34) implies a third order difference equation for the p[n] [E9|, 







1 



■p[n] {p[n-l] - p[n + l]) 



1 



-p[7i—l] (p[n— 2] — p[ra]) 



(3.42) 



which is indeed linear in p[n+l]. Substituting ( 3.40| ) gives a fifth order equation for r[ri], linear in the highest 
order term t[7i + 2]. 



3.5. Coalescence from PIV. Since the earliest works on the Painleve transcendents pi| it was known how to 
obtain the PII system from a limiting procedure or coalescence applied to the PIV, however there is in fact 
more than one such coalescence path. In fact our analysis of the scaled GUE requires the application of a second 
coalescence path rather than the one commonly employed. In the first limit the parameters (ao,ai,a2) and 
variables tiv, fflViPiVi ^iv in the PIV system scale in a way that ck2 is fixed so that 



ao 
ai 
a2 

qiv 
Piv 



V2- 



"II - ^2^" 



"11 + 1/2 , 

-e-3(a„ + 1/2) + 2^/3 



as e — > then the function qu{tii) satisfies the PII differential equation with parameter an | |l^ . 
The second limiting procedure is obtained from the first by the mapping uj introduced in (2.28) 

Proposition 20. // ai is fixed so that the variables scale like 

ao = an + 3/2 + 1/2^"^ , 
ai = -an - 1/2 , 

"2 = -y2e"' , 

tiv = e"^-2-2/3g^jj ^ 

qiV = 2^/^epu , 

2piv 
HiY + aitiv 



e-3 + 22/3e-iq„ 



(3.43) 
(3.44) 
(3.45) 
(3.46) 
(3.47) 
(3.48) 
(3.49) 



(3.50) 
(3.51) 
(3.52) 
(3.53) 
(3.54) 
(3.55) 
(3.56) 



as e — > then the function qii(tii) satisfies the PII differential equation with parameter an. Furthermore the 
third order difference equation for Hiy ( 2.8!\ 2.9C ), related to the discrete Painleve equation dPI, corresponding 
to the Backlund transformation under the shift operator Ti transforms into the third order difference equation 
for Hii ( 3.311 ), "related to the alternate discrete Painleve equation a-dPI, under this scaling. 
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Proof. Under the mapping Hiy i-^ — -ffiv, '^Piv QiVi Q^i ^ —0:2 and tjv ^ —tiv- The only equation which 
isn't immediate from the mapping is ( 3.56| ). Substituting (B.5S) for ijv and ignoring the term proportional to 
e til shows this equation is equivalent to 



Hiy = e-3(aii + 1/2) - 22/3e-ii/„, 



(3.57) 



which is precisely what results from applying the mapping lo to (3.49). The scaling of the third order difference 
equation for Hiy associated with the discrete Painleve equation, dPI, ( 2.89|) to ( |3.37 ) can be verified directly. 
□ 



4. Application to Finite GUE Matrices 



In this section we will show that the determinants in ( p.51| ) and ( 2.52 ) occur in the calculation of the quantities 
En and -Fjv, introduced in the Introduction, relating to GUE random matrices. 



4.I. Calculation of En(0; (s,oo)) and i?jv(s;a). Consider first the probability £^{0; (s,oo)). 

Proposition 21. The gap probability E]y{0; {s,oo)) is identical to the N*^ r-function of the sequence generated 
by T^^ from the corner of the Weyl chamber {ao,ai,a2) = (1,0,0), 



where the normalization of T3 [A'^] must be such that 

lim r3[iV](s;0) = 1. 



(4.1) 



(4.2) 



The resolvent kernel function R — RN(t) occurring in (1.5) is the N Hamiltonian associated with this sequence, 



RN{t)=H{t) 



(ao,ai,Q2) = (l+A',0,-JV) 



(4.3) 



Proof. From the meaning of Ejq we see from (1.1) with g{x) — e ^ that 

N 



En{0; (s, 00)) = ^ / dxi 



Introducing the Vandermonde determinant 



dxjv ]^e Y\_ i^k-Xj] 

' i=l l<j<k<N 



i<j<k<N 



(4.4) 



(4.5) 



standard manipulations following Heine (see ||3^ p. 27) allow this to be rewritten as the N x N determinant 



ENiO;{s,oo)) = ^det 



C 



■ det 



3 i j,k=0,...,N-l 

(s - xy+''e-'=' dx 



j,k=0,...,N-l 



(4.6) 
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where the second equahty follows by repeating the procedure which led to the first equality but starting from 
the Vandermonde determinant formula with Xj Xj — s (j = 1, . . . ,N). Comparing with ( ^.51 ) gives (4.1). 



RecalUng ( |2.30| ) it follows from (^) that 

En{0; (s,oo)) = exp 



(ao,ai,a2) = {l+Nfi,-N) 



dt] 



(4.7) 



and comparing with (1.5) gives (4.3). □ 

According to (2.18) R satisfies Eq. (pT^), thus rederiving the result of Tracy and Widom ( |l.5| ). We remark 
that the boundary condition follows from the fact that, with pit) denoting the eigenvalue density, 

2-N 



Pit) 



~ HN+i{t)HN-i{t) 



Considering next En{s; a) we have a generalisation of the previous case. 



(4.8) 



Proposition 22. The average EN{s;a) is the N*^ member of the T-function sequence {t3[N]{s; —a)} generated 
by T^^ from a general position on the a2 — wall of the Weyl chamber, namely (ao, «!, ct2) = (1 + a, ^a, 0), 
and so 



EN{s;a) = T3[N]{s;~a). 
The logarithmic derivative UN{t;a) is identical to the Hamiltonian associated with this sequence, 

UN{t;a) = H{t) 

{ao,ai,a2) = (l+a+N ,-a,-N) 

Proof. According to ( |l.7| ) En{s', a) is given in terms of a multiple integral by 

1 ^1 2 -I— r 

EN{s;a)^— dxi--- dxAr e^^^ (s - Xj)" Y\_ (xk-Xj)'^, 



l<j<k<N 



where the normalization C is such that 

i^Ar(s;a) - s^°(l + O(l/s))i?jv(s;0) as s 



oo. 



The method of derivation of (4.6) shows that 



1 



£'Ar(s;a) = ^ det 



(s - 2:)"+^'+'''e-^ dx 
Recalling ( |2.51| ) we thus have ( |4.9| ). It follows from this that 
EN{s;a) ^ EN{so;a)exp ( H{t) 

J So 

and consequently ( 4.10| ). □ 



j,k=0,...,N-l 



{aQ.ai,a2) = il+a+N,-a,-N) 



dt 



(4.9) 



(4.10) 



(4.11) 



(4.12) 



(4.13) 



(4.14) 



According to (2.18) UN{t;a) satisfies the nonlinear equation 

(C/;^)2 - 4(t[/^ - Un)^ + Wj,{U'^ - 2a){U'^ + 2N) = 0, 



which considering ( 4.12 ) and ( [l.5| ) is to be solved subject to the boundary condition 

UN{t;a) ^ ^(l + 0(l/<)) +i?(t). 



(4.15) 



(4.16) 
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In light of (4.8) the term R{t) decays as a Gaussian and so is negUgible with respect to the inverse powers in 
(4.16) (here we assume a ^ 0). This means the boundary condition (4.16) cannot be distinguished from the 
boundary condition with R.{t) removed; however we will see below that this latter boundary condition relates to 
a solution of ( 4.15 ) distinct from the one required in ( 4.14 ). To overcome this ambiguity we specify the t — > — oo 
behaviour of U(t) rather than the i — > oo behaviour. Now, replacing s by —s in (4.11) and changing variables 
X X + s shows 

N 



J^n[—s, a) — — e s 



/ dxi--- 
Jo Jo 



1 



1 



In combination with (4.14) this implies 

UN{t;a) 



dxi 

o(i/,s^: 

- -2Nt - 



'2xj 



l<j<*:<A' 



N{a + N) 



(4.17) 



(4.18) 



Under the action of the Backlund transformations UN(t; a) will satisfy two recurrence relations corresponding 
to the shift operators Tg"^ and Ti. For the Tg"^ sequence we relate the Painleve and GUE parameters by n = iV 
and ai — —a, and with ( p. 81 ) and (2.80) we have a third order recurrence relation in N with a fixed (suppressing 
the additional variable dependence in U) 

- N{Un+i - Un){Un - Un-i) (4t + Un+2 + Un+i - Un - C/jv-i) ^ 
{N + 1){Un+2 - Un+i) ~ N{Un " Un^i) 

- i/2(2i + Un+2 - UN)i-2a + [Un+i - C/iv)(2t + Un+i - Un))] (4.19) 
{N + 1){Un+2 - Un+i) - NIUn - Un-i) 

- i/2(2i + Un+1 - C/N-i)(2a + (Un+i ~ C/iv)(2t + Un+2 ~ C/jv-i))] = . 

In contrast the T\ sequence has the parameter correspondence n = a and a2 — —N which is an interchange of 
«! and a2 with respect to the T3 sequence. A third order recurrence relation in a, with fixed N, can be derived 
from ( 4.19 ) using the mapping ( 2.43 ) in which N t 1-^ it and then UNiit;a) 1— s- ~iUN{t; a). Alternatively 
this difference equation can be found directly from (2.89) and ( 2.90| ), with the identification (4.10). 

As noted in the Introduction, for (A^+1) x (A^+1) dimensional GUE matrices, the distribution of the largest 
eigenvalue Pniax(s) is proportional to e^* £^Ar(s;2). Hence, using (4.14), 



Pmax ( ) 



Ni^N+1 



(so) 



N^N+1 



exp I 



^2t + UNit;2)] dt 



(4.20) 



Comparison with (|l.2|), after putting 1-^ iV + 1 therein and substituting (1.5), gives an identity between 



the Hamiltonians R{t) 
verified. 



N^N+l 



and U{t; 2). Using the theory developed here this relation can be independently 



Proposition 23. The logarithmic derivative oj the average EN(t;2) is related to the resolvent kernel i?Ar(t) by 
the identity 

R'N+iit) 



UN{t]2)^2t + RN+iit) + 



RN+l{t) 



(4.21) 
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Proof. From the identifications made above we have 

UNit;2) = H{t) 
RN+i{t) = H{t) 



(JV+3,-2,-JV) 



{N+2,0,~N-1) 

and we seek to express the first in terms of the second. We note that 



H{t) 



{N+3-2,-N) 



2t + H 



{N+2,0,-N-l) 

2{N + 2) 2(iV+l) 



/o 



(4.22) 
(4.23) 

(4.24) 
(4.25) 



/2 



2{N + 2) 

To 



where the right-hand side is evaluated at {N + 2,0,— N — 1). We recognise in this expression the factors of 
H\ 



\(N+2,o-N-i) = /i(y2/o/2 - - 1) and i?'|(jv+2,o,-JV-i) = /1/2 so that it can be simphfied to 

H' 



Hit) 



= 2t + H + 

{N+3,-2,-N) H 



(4.26) 



{N+2,Q,-N~1) 



The result then follows upon making the appropriate identifications. □ 

One can verify that this transformation TiT.^^ (and its inverse) is the only nontrivial one which can map the 
Hamiltonian H into a rational function of H and H' . We have 



H 



(ao + l,Qi-2,a2 + l) 



„ (l-ai)il' 

= 2t + H+ , 
{00,01,02) rl + Cki 72 



(4.27) 



and this is rational if ai = 1 (trivial case) or if ai = (this case). All other transformations are algebraic 
functions of H and H' . 



4-.2. Calculation of Fiq{s;a). Turning our attention to Fjq, we can make the following identifications. 

Proposition 24. The average of the power of the characteristic polynomial is given by the N*^ member of the 
T-function sequence generated by the shift operator Tg"^ from the initial parameters (1 + a, —a, 0), 



FN{\;a)^UN]{\;~a), 
with the normalization of fg^-* chosen so that 

FN{\;a) - as \ ^ 00 , 

(Eq. ^.1^ )). The logarithmic derivative of the average is related to the Hamiltonian by 

VN{t;a) = H 



(l+a+N,-a,-N) 



(4.28) 



(4.29) 



(4.30) 



(Eq. Uld^)). 
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N 



Proof. We see from (1.12) that 

Fw(A;a) = -/ dxi ■ ■ ■ dx^ l[e-^Hx ~ x,r H i^k^^jf' (4-31) 

j = l l<j<k<N 

where the normaUzation is such that ( 4.29| ) is satisfied. Proceeding as in the derivation of ( |4.6D and ( 4.13| ) we 
see that this can be written in terms of determinants according to 



Fat (A; a) = ^ det 



(A - x)"+J'+'=e-^' dx 



J j,k=Q,...,N-l 

This is precisely the determinant occurring in ( 2.52| ) so we have (4.28). Recalling (2.30), we thus have 



Ar(A; a) = Fn{Xo; a) exp / VN{t; a) dt^ 



(4.32) 



(4.33) 



where VAr(i; a) is given in terms of H with [ao, ai, a2) equal to {1 + N + a, —a, —N) as in the formula (4.1C) 
relating UN{t;a) to H. □ 

Thus V^r(i;a), like UNit;a) (recall ( 4.15[ )), satisfies the nonlinear equation 



{V;^f - 4(tl4 - Vn^ + 4V^j^(y^ - 2a)iV;, + 2N) = 0. 



(4.34) 



The asymptotic behaviour ( 4.29| ) together with (4.33) implies (4.34) is to be solved subject to the boundary 
condition 



VN{t;a) ^ ^(l + 0{l/t)^ as t 



oo. 



(4.35) 



Apart from the quantity R(t), which we know decays as a Gaussian, this boundary condition is the same as ( 1.16| ). 
Thus UN{t; a) and ¥^{1; a) satisfy the same differential equation, and up to a term which decays as a Gaussian, 
the same boundary condition as t —> oo. However the t — > — oo behaviours are very different: for UN(t]a) it 
is given by ( 4.1S| ), while for VN{t;a) it is (up to a possible phase) again given by ( 4.35| ). In addition V/v(t;a) 
satisfies the AT— difference equation ( [4.19 ) and its a— difference analogue but with the appropriate boundary 
conditions. 

It was noted in the Introduction that for (A^+1) x (A^+1) dimensional GUE matrices the density, p{X) 
say, is proportional to e^'^^ F]\[{X; 2). Hence, analogous to (4.20), we see from (4.33) that 



JVi— Af+l 



On the other hand, we know p(X) 



p(Ao) 



exp 



-2t + VN{t;2)] dt 



(4.36) 



is proportional to the 2x2 determinants of (1.15). In fact FjsjlX^a), 
for general a £ Z>o can be written as an a x a determinant, a fact which can be understood in the present 
setting by considering the r- function sequence ( ^.64 ). 

Proposition 25. The average of the powers of the characteristic polynomial F]^{X;a) obey the duality relation 

FN{X;a) Fa{iX;N) 



FNiXo;a) FaitXo;N) 



(4.37) 



for alla,N e Z. 
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Proof. First, note that by reversing the steps which led to (4.6) and recalling (4.31) the determinant of (2.64) 
which specifies fi [n] can be written as a multiple integral to give 

TiM(A;a2) = CFn{iX;~a2) , (4.38) 

where |C| = 1. On the other hand, we see from ( ]2.13| ) that the Hamiltonian corresponding to Ti[n], 



H 



(ao,cx.i,a2) — {l — <y-2-\-' 



n,—n,OL- 



2) 



satisfies the differential equation ( 4.15 ) with 



N 



-OL2, 



(4.39) 



(4.40) 



(Note that for the latter identification to be possible, we require a £ Z>o.) Furthermore, from (4.38) and (4.2£), 
we have that 



(4.41) 



which with the substitutions ( |4.40| ) is identical to the boundary condition of ( |4.35| ). It follows from these facts 
that 

Fa{i\ N) = Fa{i\^; N) exp ( j VN{t; a) di) . (4.42) 



Comparison with (4.33) then yields (4.37). □ 

From the definition (4.31) this identity implies the integral identity 



/oo i^oo ^* 

dxi--- dxAT Jl e"''?(A - Xj)° Y[ [xk-Xjf 

j=l l<j<k<N 

/OO poo Q 

dxi--- dxa ]^e"^?(A- ixj)^ J| {xk ~ xj)"^ 
-oo J —CO i^.-^i,^„ 



and from (4.32) gives the determinant identity 



det 



(A - xYx^+'^e-''^ dx 



j,k=0,...,N-l 



3 = i 



Cdet 



l<j<k<a 



(A - i2;)^x^+'^e-^' dx 



j,k=0,...,a~l 



(4.43) 



(4.44) 



The integral identity ( 4.43 ) has been derived earlier in the context of a theory of generalized Hermite polynomials 
based on symmetric Jack polynomials Q , and in fact can be generalized so that on the left-hand side the exponent 
2 in the product of differences is replaced by 2c and i-^ ca;| in the Gaussian, while on the right-hand side 
this same exponent is replaced by 2/c. Regarding the determinant identity, noting that the right-hand side is 
proportional to 



det 



(A-ix)^+^'+'=e-"'dx' 



HN+j+kW 



j,k=0,...,a-l 



(4.45) 



= Cdet 

j,fc=0,...,a-l 

this gives a determinant formula for Fn{X; a), equivalent to that given by Brezin and Hikami Q. Determinants 
with a Hankel structure constructed with orthogonal polynomial elements are termed Turanians and their 
positivity and other properties such as relations with novel Wronskians have been investigated extensively, and 
reviewed in Karlin and Szego . Explicit evaluations of Turanians of the Hermite polynomials in terms of the 
Barnes G- function where the initial degree of the polynomial is zero {N — 0) have been given by Radoux . 
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4-3. UN{t]a) and VN{t;a) as Painleve transcendents. The formula (4.10) relates U]\[(t;a) to the Hamiltonian 
( 2.13 ). Substituting the appropriate values of ai and 0:2 in the first equality of ( ^.13 ) and recalling ( 2.16| ) shows 



UN{t;a) 



{q'f + 2t){q^ + 2tq - 4a) - ^ + Nq. 



8 



2q 



where q satisfies the PIV equation (2.1) with 



a = 2^^ + 1 



f3- 



(4.46) 



(4.47) 



In the case a = the functional expression ( 4.46| ) agrees with that presented earlier [^|3^, although the 
transcendent q in the earlier works is the PIV transcendent with a — 2N — 1, f3 — rather than a — 2N + 1, 
/? = as given by ( 4.47 ). In fact it follows from the work Q that in general Eq. ( |2.18| ) has more than one 
expression in terms of Painleve transcendents. In the case ai = iV, Q!2 = of this equation the results of M give 



the functional expression implied ( 4.46|) with q the PIV transcendent specified with the parameters /? = and 
either a = 2N + 1 or a = 2N — 1 thus reconciling (4.46) in the case a = with the results of We remark 

that the theory of H gives distinct functional forms for the derivative, 



U'r,it;0)^-^eq'~lq' 



tq, 



(4.48) 



in the two case a = 2N + e (e — ±1) 



Comparison of the formulas (4.30) and (4.10) shows V/v(t;a) is given by the same Hamiltonian as UN{t',a). 



Thus (4.46) remains true with the function Unit] a) replaced by V/v(i; a) on the left-hand side. 



5. Edge Scaling in the GUE 

5.1. Calculation oJE^°^^{s) and E^°^^ [s\ a) . To leading order the support of the eigenvalue density for x GUE 
matrices is the interval {—^/2N, y/2N). To study distributions in the neighbourhood of the largest eigenvalue 
one shifts the origin to the edge at V2iV and then scales the coordinate so as to make the spacings of order 
unity in the N —> 00 limit. This is achieved by the mapping [p^ 

A 



Ai 



'2N - 



v^A^ 1/6 ■ 



(5.1) 



Suppose we make this replacement (in the s-variable) in the probability E]y{0; (s, 00)) as specified by (4/7). Then 
with 



£;"°"(s) := lim En(0;V2N 



(5.2) 



\/2iVi/6. 

(because the eigenvalue density is not strictly zero outside the leading order of its support the edge is referred 
to as a soft edge) we see that 



E' 



^f*(s)=exp(- r r{t)dt) , 



(5.3) 



where 



r(t) — lim — ^ 



R{V2N 



^/2N^/i 



(5.4) 
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Furthermore, it follows from changing variables s t-^ \/2N + s/V2N^^^ in ( |l.3| ), replacing R{\/2N + s / ^/2N'^/^) 
by \/2N^/^r{s) and taking the limit N ^ oo that r{s) satisfies the differential equation 



(r")^+V((r')^-sr' + r) = 0, 



(5.5) 



a result first obtained by Tracy and Widom [Q. Equation (5.5) is a particular case of the Jimbo-Miwa-Okamoto 
a form of the Painleve II equation. We will find that the edge scaling is essentially the second coalescence limit 
of the PIV system to the PII as discussed in Subsect. 3.5. 



Proposition 26. Define the scaling limit of the quantity Epj{s;a) by 



E'°^\s;a) := lim (Ce-'"''^^ENii 



Then 



where 



E'°f\s; a) = E'°^\so; a) exp ( / u{t; a) dt^ 



u(t:a)= lim —( — at + UN(t:a)] 



V2N 

.2l/3^(_2l/3i) 



/2lV+t/y2Afi/6 



{ao,ai) = {l-a,a) 



(5.6) 

(5.7) 

(5.8) 
(5.9) 



with H{t) is given by (3.4). The function u{s;a) satisfies a second order second degree differential equation of 
the general Jimbo-Miwa- Okamoto a form of the Painleve II equation 



{u"f + Au'({u'f - su' + - a2 = 0, 



subject to the boundary condition 

u{s; a) 



(5.10) 



(5.11) 



The function u{s; a) also satisfies the third order difference equation, related to the alternate discrete Painleve 
a-dPI equation, equation 

a a+l 



~ s — [u{s; a+l) — u{s; a)] 



(5.12) 



u(s; a+l) — u(s; a — 1) u{s; a+2) ~ u{s; a) 

Proof. Unlike the probability En{0; {s, oo)), we do not expect the soft edge scaling limit of the quantities En{s; a) 
as specified by (4.11) to be well defined. For example, in the case a = 2, it is the combination e^^^ En{s;2) 
which is proportional to Pmax(s), and thus which should have a well defined scaling limit. This suggests that for 
general a we consider the scaling limit of 



Ce-'^^'/2^Ar(s;a). 



According to (4.14) we have 

e-'^«V2^^(g. ^ e-''"o/2£;jv(so; a) exp ( ( [-at + UN{t\ a)] dt^ . 



(5.13) 



(5.14) 
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Relation (5.8) follows from the definitions of £'*'°^*(s;a) and u{t;a) and (5.14). The scaling in ( pT^ ) is identical 
to the coalescence of the PIV system to the PII defined in Prop. 20, with the identifications e = {2N)~^/^, 
ct(iv)i = and ( 4.10| ) for the relationship of the PIV Hamiltonian and UN{t;a), and the scale changes 
til = ~2^^^t and Hu{tu) — —2^^^^u{t;a). Proceeding as in the derivation of (5.5) we find from 



(ao,Qi) = (l-a,a) 



( 4.15 ), or from ( 3.17 ) using the change of scale ( |5.9| ), that u satisfies ( 5.10 ). 

To formulate the boundary condition for u{t; a), we first recall |Q that the s 
r(s) in (5.5) is given by 



r(s) 



s— > — oo 



2 1 

V4S 



6As* 



O 



(^). 



-00 boundary condition of 



(5.15) 



and this corresponds to the asymptotic behaviour 



E^°'\s) ^ exp(i- + llog(- 
Also, we know that -E'soft('S; 2) is proportional to the derivative of E^°^^{s), which implies 

o3 



This suggests that for general a we seek a solution of ( ^.10 ) with the s —^ —00 boundary condition 



u(s; a) 



i=i 



Substitution shows that in fact (5.1C) has a unique solution of this form, and furthermore 



4a2-l (4a2-i)(4a2-9) 
8 64 



(5.16) 



(5.17) 



(5.18) 



(5.19) 



The third order difference equation is just ( 3.37 ) with the scale change ( p.9[ ). □ 

We see from (5.11) that the asymptotic expansion of u(s;a) terminates for a equal to half an odd integer. 
Recalling (5.!;) this is the case of ai half an odd integer of the PII Hamiltonian (3.4). From the text about (3.29) 
we know that this is precisely the parameter value for which the PII Hamiltonian can be expressed as a rational 
function of t. 

The fact that E]y{s; 2) is proportional to Pmax(s) implies that the corresponding quantities in the scaled 
limit, E^°^^{s; 2) and PmaxC'*) are proportional. Hence 



Pmax(s) = Pmax(so) exp ( T u{t; 2) dt) . 



The relation 



Pmax('*) 



E''°^\s) 



ds 



(5.20) 



(5.21) 



(Eq. ( |l.2|) ) with E^°^^{s) specified by (5^) then implies an identity between transcendents analogous to that of 
Prop. |2|. 



Proposition 27. The quantity u(t;2) of (5.2C) and the quantity r(t) of (5.5) are related by 

uit;2)^ j^\ogr{t)+r{t). 



(5.22) 
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Proof. With 



H{t) 



H[n], 



it follows from (5_^) and the substitution r{t) = u{t;0) that ( ^.22| ) is equivalent to 



H[2]^-\ogH[0]+H[0]. 



But from the analogue of the first equality in (3.21), together with (3.15) and Table |[ 



H[2] = T^H[0] = H[0] + q[0] + + T2 [ q[0] + 



H[0] 



/i[0] 
1 



ai 

m 



cti=0 



Furthermore the first equality in (3.18) together with (3.4) in the case ai = give 



2g2[0] +t-p[0] dt 



= -logH[0]. 



Note. According to (5.9) we therefore have 

-2-i/3,.(-2-i/3<;a) = 1/2 [q'{t,a-y,)]^ - 



g2(t,a-y2) + ^ 



- aq{t,a-i/2), 



(5.23) 



(5.24) 



(5.25) 



(5.26) 



where q — q{t,a) satisfies the PII equation (3J). Also the first member of ( |3.18|) , along with (3^) and (3.€) 
gives 



-2-l/3|,(_2-V3,.„)^_y^ 



q'{t,a-i/2)+qHt,a-i/2) + - 



In the special case a = we see from (3.7) that ( [3.27 ) simplifies to read 

u'it;0)^~q^it,0) 



while (5.26) simplifies to read 

u{t;0) = [q'{t,0)f -tq\t,0)-q\t,0). 
The results ( ^^ and (|^, deduced in a different way, can be found in |3 



(5.27) 

(5.28) 
(5.29) 



5.2. Calculation of F {X;a). Let us next consider the scaled quantity 



F^°f*(A;a):= lim (Ce-''^'/^FN{X;a)) 

N—>-OG \ / 



_ , (5.30) 

where Fjv(A;a) is specified by (4.31). Because of the analogy with iJ'^°^*(s; a), which follows from the identical 



structure of ( [4.141 ), ( |4.15| ) and ( [4.33| ), ( |4.34| ), analogous to ( |5.7D we have 

F''°f*(A; a) := F'°^\\o; a) exp ( / v{t; a) dt) , 



where 



v{t; a) ~ lim 



at + V/v(i; a) 



iM\/2]V+t/\/27Vi/6 



(5.31) 



(5.32) 



satisfies the differential equation ( 5.10| ) and the difference equation ( 5.12 ). The only difference between the 
logarithmic derivatives v(t\a) and u{t;a) is the boundary condition. 
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Proposition 28. The scaled averages of the powers of the characteristic polynomial F^°^^{X; a) for a G Z>o have 
the determinantal form 



F''°f*(A;a) = {-1)"^''-^^/^ det 



r 

dXP^ 



Ai(A) 



j,k=Q,...,a-l 



(5.33) 



which was shown by Okamoto ^26^ to he a r-function sequence for Painleve II (recall Prop. 18). Furthermore 
this scaled average has a multiple integral representation of the Kontsevich form JT^ 

(_l)a(a-l)/2 



dvi 



The logarithmic derivative v(t; a) has the asymptotic expansion 



dv,\{e</^-^^^ n K-^.f 

j=l l<j<k<a 



v(t; a) 



t-oo it ■ 32t5/2 



(5.34) 



(5.35) 



Proof. For positive integer a we can determine the A — > oo behaviour of F^°^^{X; a), and thus the corresponding 
behaviour of v{t;a), by making use of the scaled form of the right-hand side of ( [4.31 ). To determine this scaled 
form, we first require the explicit values of the constants in ( 4.31 ), ( 4.44 ) and (^.30|) . Let us denote these constants 
by Ci, C2, C3 respectively. Then from (4.31) and (4.29) we read off that 

N 

dxM JJ^ e~ 

2-JvV2(27r)^/2G(Ar+2), 



Ci :=ciN) 



dx^ ■ ■ 



l<j<*:<A' 

(5.36) 

where G{x) denotes the Barnes G-function, characterised for x a positive integer by the functional property 
G{x + 1) = r{x)G{x) and the initial value G(l) = 1. The integral evaluation in ( 5.3(: ) can be derived by making 
use of the Vandermonde determinant identity ( [4.5| ) written in terms of Hermite polynomials. 

The proportionality constant C in ( 4.44 ) is the same as that in ( 4.45 ) and thus from ( 5.36 ) given by 

c{N) 



C2 = 



c(a) 



(5.37) 



Finally, we seek the value of the constant in (5.30). We know that in the case a = 2, e "'^ ^'^FN{X;a) is 
proportional to the eigenvalue density for (A^+1) x (A^ + 1) dimensional GUE matrices. Specifically 



P(A) 



(N + l) 



c{N) 



-e^^ Fn{X;2). 



(5.38) 



N^N+l ~ 'c{N+l) 

Since it is the combination p{X) dX which has a scaled limit, it follows that in the case a = 2, C3 = (A^ + 



l)c{N) / {c{N + 1)V2N^/^). This suggests that for general a e Z>o we should choose 



G3 = Ca 



r{N+l + a/2) c{N) 



1 



F{N + 1) c{N+a/2)\^N^/e, 

where p is a power to be determined and Cq depends only on a (p = 1, Cq = 1 for a = 2). 

Substituting (|36|), dS^Tl ), ( ^39| ) in ( p3l| ), ( ^li] ), ( |53o| ) respectively we see that for a e Z>o, 

FiN + l+a/2) c{N) / 1 

-a/2)c(a)\J2m/<^^ 



(5.39) 



F"°"(A;a) =Ca lim 



X det 



F{N + 1) c{N + a/2)c{a)\y/2m 

(A - ix)^xi+^e^''^ dx 

J Ah^V2lV+A/\/2Afi/6 



i,fc=0,...,a-l 



(5.40) 
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But analogous to the equality in (pq) we have 



det 



j,fe=0,...,a-l 



= det [ / (A - ^xf+^+'e-^" dx 



j,k=0,...,a~l 

(5.41) 



This can be further rewritten by noting that analogous to ( p. 60 ), 



(A - ixj^+^+'^e-^' dx = (-2)-(^+'=)e^'^-^(e-^' / (A - ixfe'^" dx) 



= (-2)-(^-+'^)e^^2--V^-;^(e-^F.(A)). 
Making use of the asymptotic expansion for the Barnes G-function 

logG(x + l) - ^loga;-^x2 + ^log2^--^loga; + 0(l), 
and the Plancherel-Rotach asymptotic expansion of the Hermite polynomials 

exp{~xy2)HN{x) = ^-3/42A'/2+l/4(^,)l/2^-l/12|^Ai(_^/3l/3) _^ 0(iV-2/3)| 



(5.42) 



(5.43) 



(5.44) 



where x = {2N)^/'^ — 2^^/^3^^/^A^~^/^i and with Ai(a;) denoting the Airy function, we see from Eqs. (5.4C), 
(5.41) and ( [S.42| ) that with p — a/2 in (5.39) and appropriate Ca, the determinantal representation ( |3.33[ ) holds. 
Furthermore, in the case a = 2 we read off the functional form 



(Ai'(x))'- Ai(2:)Ai"(x) 



(5.45) 



which is the known expression for the scaled soft edge density in the GUE. Another point of interest, which 
follows from the integral formula 



Mix) f 



exp (y^'v^ ~ XV ) dv, 



is that ( 5.33| ) can be written 



F^°ft(A; a) = ^ ' ^ det 

^ ^ (-27ri)'' 



exp 



V3 



3 - Xv]v^+''dv 



j,fc=0,...,a-l 



(5.46) 



(5.47) 



Thus, reversing the reasoning leading from (4^) to (4^) we have the multiple integral representation ( 5.34 ) for 
F^°^^{X;a), which is an example of the class of integrals studied by Kontsevich |]l8| . 

Consider now the asymptotic form of (5.33). In the case a = 1 this is just the Airy function, which has the 
known x ^ oo asymptotic form (see e.g. po| p. 116) 



Ai (x 



fc=o ^ 



where ^ := Uq — 1 and 



(2fc+l)(2/c + 3)---(6/c-l) , ^ ^ 

Ul. = ■ --; : , K > 1. 

(216)'=/c! 



(5.48) 



(5.49) 



Application of the T-Function Theory of PIV, PII to the GUE 



37 



It follows from this and ( ^.33 ) that for general a G Z>o 

2a 



logf^-^A;a) ^ ^ -^A3/2 + Clog A + co + ^ ^ 



A— >oo 3 



5z_ 

X3J/2 



(5.50) 



which in combination with ( ^.31 ) implies that we must seek a solution of ( 5.10| ) (with u replaced by v) subject 
to the boundary condition 



.(t;a)^^^-a^^/^ + f + ;^^. 



(5.51) 



Substitution of ( 5.51 ) in ( 5.10 ) shows there is a unique solution of this form, with 

a2 a(l+4a2) 



(5.52) 



given by ( ^.35 ). □ 



6. Conclusions - A Programme 

We have applied the Okamoto r-function theory of PIV and PII to the computation of En^s^o) and FN{s;a) 
for the GUE and its scaled soft edge limit. As noted in the Introduction, the Okamoto r-function theory applies 
equally as well to the computation of Em{s; a) and Fn{s:, a) for all matrix ensembles with a unitary symmetry 
and classical weight functions (p^). Thus we expect to be able to compute £'Ar(s; a) and F]^{s; a) in the cases of 
the Laguerre, Jacobi and Cauchy ensembles (special cases of Fisf{s; a) have been evaluated in terms of Painleve 
transcendents for the Laguerre ensemble , and for the Jacobi ensemble j|] ) . In future studies we will undertake 
this task by following the programme used here for the GUE, the main steps of which can be itemised as follows: 

— From the definitions of the gap probability En{0',I), I a single interval including the boundary of the eigen- 
value support, and En{s; a), Fn{s; a) as A-dimensional multidimensional integrals they can be converted into 



N X N determinants analogous to (4.6), (4.13) and (4.32) respectively. 



— Using an identity analogous to (2.60), the determinants can be put into the double Wronskian form ( ^.58| ), 
with d/dt replaced by 

d , , d , , 

'df '^'-'^dt (^-1) 

in the Laguerre and Jacobi ensembles respectively. 

— The Okamoto r-function theory of PV and PVI [^|2^ gives these same determinants as r-function sequences, 
in which the initial members are r[0] = 1, and r[l] the solution of the particular classical equation associated 
with relevant Painleve transcendent when the parameter sequences begin on a wall of the Weyl chamber in 
the affine space of parameters. The classical solutions, and their polynomial specialisations, are noted for each 
of the Painleve transcendents in Table ra, 
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PJ 


Classical Solution 


Classical 
Orthogonal Polynomial 


PI 






PII 


Airy 




PHI 


Bessel 




PIV 


Hermite- Weber 


Hermite 


PV 


Confluent Hypergeometric 


Laguerre 


PVI 


GauB Hypergeometric 


Jacobi 



Table 8. Classical solutions of the Painleve transcendents. 



— The logarithmic derivatives (with d/dt replaced by (|6.l| ) as appropriate) R]y{s),UN{s),VN{s) coincide with 
the Hamiltonians in the Painleve theory and as such satisfy certain second order second degree ODEs of the 
Painleve type. 

— The T-function sequence {ro[A^](t; a)}jv>o, say corresponding to Fpf{s;a), is simply related to another r- 
function sequence {Ti[a]{t; N)}a>o- Both r- functions relate to the same Hamiltonian but result from the 
action of different shift operators. Because the shifts are commutative one has 

ro[N]{t;a) ^ T,[a]{t- N) 

To[N]{to;a) Ma]{to;N) ' ^ ' ^ 

Identities of this type for the Laguerre and Jacobi ensembles, written as multiple integrals, are already known 
from 0. 

— For all the independent shift operators and sequences of (7[n],p[n], iJ[n], r[n] there exist difference equations 
generated by the Backlund transformations of these shifts. It has been conjectured that all the difference 
equations arising in this way are discrete Painleve equations satisfying integrability criteria such as singularity 
confinement analogous to the Painleve criteria. 

— In the appropriate edge scaling limit, the analogues of r{s),u{s),v{s) are Hamiltonian functions for PII or 
PHI, and satisfy the corresponding second order second degree equation. 
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